Draft version February 5, 2008 

Preprint typeset using I4TgX style emulateapj v. 6/22/04 



STAR FORMATION IN DISTANT RED GALAXIES: SPITZER OBSERVATIONS IN THE HUBBLE DEEP FIELD SOUTH 
Tracy M.A. Webb', Pieter van Dokkum^'^, Eiichi Egami Giovanni Fazio^, Marijn Franx', Eric Gawiser'-'-^, David 

HERRERA-•^ JIASHENG HUANG^, IVO LABBE^, PAULINA LIRA^, DANILO MARCHESIN^•^ JOSE MAZa'', RYAN QUADRI^, GREGORY 

Rudnick', AND Paul van der Were' 

Draft version February 5, 2008 



in 
o 
o 

> 
O 

o 



> 
oo 



in 
o 

Oh; 

6 



ABSTRACT 

We present Spitzer 24/im imaging of 1.5 < z < 2.5 Distant Red Galaxies (DRGs) in the lO'x 10' Extended 
Hubble Deep Field South of the Multiwavelength Survey by Yale-Chile. We detect 65% of the DRGs with 
Kab < 23.2 mag at S24^m ^ 40/iJy, and conclude that the bulk of the DRG population are dusty active galaxies. 
A mid-infrared (MIR) color analysis with IRAC data suggests that the MIR fluxes are not dominated by buried 
AGN, and we interpret the high detection rate as evidence for a high average star formation rate of < SFR > = 
130±30 MQyr"'. From this, we infer that DRGs are important contributors to the cosmic star formation rate 
density at z 2, at a level of ~ 0.02 M0yr"'Mpc"-' to our completeness limit of Kab = 22.9 mag. 

Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: starburst — infrared: galaxies — 
(ISM:) dust, extinction 



1. INTRODUCTION 

Through a variety of observational methods we are build- 
ing a census of the high-redshift universe, and can now di- 
rectly study galaxies in the process of forming. Color selec- 
tion criteria effectively select large sampl es of galaxies at z > 
1 such as Lyman-break galaxies (LBGs) JSteidel e t alJ 19961). 
and th e massive Distant Red Galaxies (DRGs) dR^nxetall 
UdO?). Selected by rest-frame optical colors to lie at z > 2, 
DRGs exhibit properties of both passively evolved stellar pop- 
ulations and dusty starburst galaxie s ( Forster Schreiber et al. 
l2()04l iLabbe et al.ll2(M iPapovich et al...200 5). Characteriz- 
ing this complex population and placing it into the context of 
the star formation history of the universe will elucidate our 
overall understanding of the assembly of massive galaxies. 

Recent advances in the capabilities of infrared facilities 
have opened a new window onto the high-redshift universe. 
Through the direct detection of dust enshrouded activity, mid- 
and far-infrared (M/FIR) studies provide orthogonal informa- 
tion to that gathered in the ultraviolet (UV), optical and near- 
infrared (NIR), which can be heavily biased by dust extinc- 
tion. The work presented here addresses the nature of z ~ 2 
DRGs through Spitzer MIPS-24^m and IRAC imaging. We 
assume a Hm = 0.3, J^a = 0.7 cosmology and Ho = 70 km s"' 
Mpc"' throughout. 

2. THE SAMPLE, SPITZER OBSERVATIONS, AND PHOTOMETRY 

The DRG sample was drawn from the lO'xlO' Extended 
Hubble Deep Field South (EHD F-S) of the Multiwav elength 
Survey by Yale-Chile (MUS YC) JGawiser et al.l2005t) DRGs 
are defined by (J - K)yega. > 2.3 mag jFranx et all l200l 
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Ivan Dokk um et a and, in this paper, we focus on 

DRGs with Kab < 23.2 mag (total magnitudes are given 
throughout), corresponding to the MUSYC 50% complete- 
ness limit (Quadri et al., in preparation). 

Photometric redshifts were derived from UBVRIz'JHK 
photometry using the code presented in Rudnick et al. (200i1, 
120031) . through linear combinations of galaxy templates, with 
an accuracy of ( \zspec-Zphor \ / C^+Zspec) ) =0.05 forz > 1.5. 
Here, we have restricted our analysis to redshifts 1.5 < z < 
2.5 where the 6.2/im, 7.7/im, and 8.6/im Polycyclic Aromatic 
Hydrocarbon (PAH) features of star-forming galaxies fall into 
the Spitzer 24/im filter. This provides a sample of 79 DRGs 
with a median redshift of z = 2.0. 

The 24/im observations were taken in the MIPS photometry 
mode and consist of 6 separate pointings. The central deepest 
pointing represents ^^1 hour of frame time, and the 5 flanking 
pointings have 35 minutes of frame time each. The raw data 
were reduced and combined with the Data Analysis Tool de- 
veloped by the MIPS instrument team ( Gordon et al. 2004). 
The final image covers ^104 arcmin^. 

Imaging at 3.0/im, 4.5/im, 5.8/im and 8.0/im was taken in 
the IRAC mapping mode and mosaicked to produce a ^140 
arcmin^ image. Individual frames of 200 seconds were com- 
bined for a final total frame time, per location on the map, 
of 20 minutes. The Spitzer Science Center provides a Ba- 
sic Calibrated Data product with flat-field corrections, dark 
subtraction, and linearity and flux calibrations. Additional 
steps conducted with the IRAC team customized software 
included pointing refinemen t, distortion correction, and mo- 
saicking ('Huang et alJ2004l) . 

Photometry on the 24/im image was performed using 
the point spre ad function (PSF) fitting program DAOPHOT 
llStetsonll1987h . The PSF was characterized using bright iso- 
lated point sources and then used to iteratively fit and sub- 
tract each object in the image. A small fraction (10/79) DRGs 
were highly confused with neighboring objects and these were 
removed from further analysis. We reached an average rms 
depth of ~13/iJy, and fou nd the sourc e counts to be in good 
agreement with Papovich et alJ (|2003)- Because not all the 
DRGs are detected at 24/im, their flux densities and upper 
limits were determined in a consistent manner by performing 
aperture photometry at each /T-determined position after first 
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subtracting all other objects in the catalog. We used 6" di- 
ameter apertures, and used the PSF stars to apply an aperture 
correction factor to total flux densities at 40". 

Due to the smaller beam size at the shorter IRAC wave- 
lengths (< 2"), and because highly confused DRGs have been 
removed from the sample, PSF fitting was not employed for 
the IRAC measurements. Aperture photometry (3" diame- 
ter) was performed at the location of each DRG for each 
IRAC channel, again using isolated point sources to deter- 
mine aperture corrections to total fluxes at 12.2". The average 
rms depths are 2.6/iJy, 2.7/iJy, 4.3/iJy and 4.0/iJy at 3.6/im, 
4.5/im, 5.8/im, and 8.0/im, respectively. 

3. MIR CONSTRAINTS ON THE NATURE OF DRGS 

At 1.5 < z < 2.5, the 24/im filter samples ^ 6-10/im in 
the rest-frame and the 24/im flux offers a powerful method 
of differentiating between two basic SED types: active dusty 
galaxies, i.e., those that are powered by either star formation 
or active galactic nuclei (AGN), both of which produce sub- 
stantial mid-infrared emission, and passively evolved systems 
whose stellar flux drops drastically longward of ^l^im (see 
Fig.Q. In pure starburst galaxies, the MIR emission is dom- 
inated by PAH features which are strong relative to the un- 
derlying dust continuum, which only begins to rise above 
10/im. The hard radiation field of an AGN destroys PAH car- 
riers and the continuum emission from hot small dust grains 
is strong throughout the MIR (Genzel & Cersarskv 2000). 

In Fig. n we show the 24/im flux densities for the DRGs, 
normalized to their rest-frame Ly, and overlaid with three 
template SEDs; 65% are detected at 24/im and the remain- 
der are shown as 3<t upper limits. This detection rate is in 
line with that found by (Papovich et al. 2005) for DRGs the 
GOODS Chandra Deep Field South, to a similar depth. All 
of the 24/im-detected DRGs are incompatible with old stel- 
lar populations and dusty starburst or AGN galaxies are re- 
quired to produce such strong mid-IR emission. The sample 
shows a wide range in 24/im flux which we interpret primar- 
ily as a range in star formation rates. However, a real dif- 
ference in PAH strengths, relative to the continuum, due to 
a range in metallicity or the hardness of the radiation field 
(e.g., Engelbracht et al. 2005; Hogg et al. 2005) could also be 
important, as could a variation in the level of AGN contami- 
nation of the MIR flux, though as we discuss below we deem 
this latter possibility unlikely. 

The SEDs of DRGs which are not detected at 24/im (35%) 
are consistent with passively evolved systems, but also lower 
luminosity or less dusty starbursts, or starbursts which have 
weak or absent PAH signatures as explained above. To inves- 
tigate the nature of the 24/im-faint DRGs we used a stacking 
analysis to measure their mean 24/im emission and find 524^m 
= 10 ± 2 /iJy. Thus, the remaining 35% of the sample can- 
not be entirely composed of passively evolved systems but 
contains some fraction of dusty starburst or AGN galaxies. 
A bootstrapping analysis indicates that this measurement is 
not dominated by a small number of systems, but beyond this 
we cannot constrain the fraction of starburst or AGN systems 
within this group. 

Combining the 24/im measurements with IRAC da ta pro- 
vides a color diagno stic of AGN activity (Fig . 13 ( Iviso n et alJ 
I2OO4I lEgamiet alJ l2(K)l iHuang et alJ UmS). This tech- 
nique reUes on the different rest-frame NIR SED proper- 
ties of AGN and starburst galaxies: stellar emission in star- 
burst galaxies leads to a relatively flat SED in the rest-frame 
NIR, while hot dust around an AGN produces a power- 
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Fig. 1. — The 24/im fluxes or limits of the DRGs, normalized to their 
rest-frame Ly. Points are DRGs with 24/im detections, and arrows show 3cr 
upper-limits for the remainder. The single solid point below the arrows de- 
notes the stacked measurement of DRGs which are not individually detected. 
Three SEDs are shown: the green line corresponds to an elliptical galaxy 
model I Coleman, Wu, & Weedman 1980) extrapolated to longer wavelengths 
using ISO observations of local ellipticals; the blue line to the observed SED 
of Arp220; and the red line to the observed SED of NGC1()68 (Jiasheng 
Huang, personal communication). Because of the rapid variation in flux with 
wavelength for the Arp 220 SED in the MIR we also show this SED smoothed 
by the 24/im filter transmission curve (smooth black line). 



law that increses toward longer wavelengths. We have de- 
tected 50% of flie 24/tm-bright DRGs at 4.5/im and 8/im, 
another 20% at 4.5/im, but not 8.0/im, 5% in neither, and 
25% faU off the IRAC channel 2/4 field-of-view. The col- 
ors (Fig. 2} indicate that star formation dominates the MIR 
emission of the DRGs, but as a conservative cut we con- 
sider the four objets with 58^n]/5'4.5pm > 1.5 to have possi- 
ble AGN contamination. This is in agreement with results 
from the X -ray (Rubin et al. 2004; Reddy et al. 2005) and op- 
tical (|van Dokkum et al. 2003) which also indicate that AGN 
do not significantly contribute to the observed luminosity of 
DRGs. It should be remembered, however, that very heavily 
obscured or weak AGN, though energetically unimportant in 
terms of observed MIR properties, may still be present. 

4. STAR FORMATION RATES OF DRGS 

PAH and MIR emission can be used to esti mate the 
current star formation rate (SFR) of a galaxy SWu et all 
2005). Using an Arp220 SED template we extrapolated from 
24fim/(l+z) to rest-frame 6.75/im. We then estimated Lir 
through the 6.75/im-FFiR relation calibrated locally with the 
Infrared Space Ob servatory, which shows a scatter of ^50% 
("Elbaz et alJ2002 ). This value can then be converted to a SFR 
following the Lir-SFR relationship of Kennicutt ( 1998) with 
an uncertainty of a factor ^ 2-3. We note that adopting an 
M82-like SED instead of Arp220 in the first step results in 
SFR estimates which are systematically 1 .4 x lower 

In Fig. |3lwe plot the 24/im star formation rates or upper 
limits of the individual DRGs with redshift. The population is 
not homogeneous and individual galaxies range in SFR from 
< 30 M0yr"' to --1000 MqJV'K A strong AGN contribution 
to the MIR flux will contaminate this estimate, by boosting 
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Fig. 2. — A MIR color-color plot. The tracks correspond to two starburst 
SEDs (left-most, orange-Arp220; red-M82) and three AGN dominated SEDs 
(right-most, magenta-MRK23 1 ; puiple-Seyfert 1 NGC5506; green-Seyfert 2 
NGC1068). The SED tracks begin at z = 1.5 and are marked by circles mz = 
0.5 steps, to z = 3.5. The open points show DRGs from this work which are 
detected at 24/im, with the radio detected objects indicated by the large red 
circles. Also shown are Spitzer observed SM Gs (solid circles: [ Egami et alj 
120041) and optically-faint AGN (solid triangles: lRigbv et all2003) . 



the inferred SFR, and sources which may contain AGN are 
marked. Also highhghted are three DRGs wh ich are detected 
in the deep radio map of Huvnh et alJ J2005 '). two of which 
are possible AGN based on their MIR colors. Studies indicate 
that the fraction of AGN contaminati on in the uJy radio pop- 
ulation increases at higher redshift iRichardsll200Ct iGarretll 
1^01), and thus, the radio detected DRGs may indeed contain 
weak AGN. However, we note that these systems fall within 
1(7 of the FIR/radio flux correlation of Condon ( 1992) (with 
no systematic offset), which only holds for starburst galaxies, 
and thus significant AGN contamination is unlikely. 

Using a stacking analysis, we calculated an average SFR 
for the entire DRG sample of <SFR> = 130±30 Mgyr"'. 
This is the uncertainty-weighted mean of all 69 DRGs, as- 
suming no AGN contamination and a mean redshift of z = 2, 
with the uncertainty in the mean estimated from a bootstrap- 
ping analysis. This result is in good agreement with the opti- 
cal SED modeling results of Forster Sch reiber et alJ ( |2004» of 
120 M0yr"', and X-ray and submillimeter stacking analyses 
which estimate ^100-200 M^^^vr"' and ^ 100 Mm vr"', respec- 
tively dRubinet al...2004i iReddv et all 12005; .Knudsen et all 

A significant source of error in this calculation is the pho- 
tometric redshift uncertainty, since the rest-frame flux of star- 
burst galaxies changes rapidly with wavelength in the region 
surrounding the PAH features. This effect was quantified by 
recalculating SFRs for simulated redshifts within Gaussian 
confidence intervals, and the resulting uncertainties are shown 
in Fig. |3l The uncertainties are substantial, ranging from a 
factor of two to almost an order of magnitude. 

We calculated the optically normalized star formation rates 
(SFR/Lv) for the DRGs, which provide an indication of their 



Fig. 3. — The SFRs for DRGs as a function of redshift (circles). Filled cir- 
cles mark DRGs which may contain AGN. and red circles mark radio detected 
galaxies. The solid line con'esponds to the average effective SFR depth as a 
function of redshift, for these data. Also shown for reference is the median 
SFR for radio-detected submm galaxies from|Chagman e^^^ 1 . 200 5) (solid 
star) and that of Lyman-break galaxies from lShap!evT^^l200lF foDen star). 



age and/or dust content. We find < SFR/Lv 19 , which is 
in excellent agreement with the submillimeter study of DRGs 
by Knudsen etal. (2005), who measure < SFR/Ly >^ 20. 
However, the values for the DRGs show large scatter, ranging 
from SFR/Ly < 5 to 200. Monte Carlo simulations indicate 
that the redshift uncertainties cannot account for the full range 
of this scatter and, moreover, the three radio detected galaxies 
have the highest SFR/Ly, for a given Ly; as their radio flux 
confirms that they are extreme systems, their high SFR/Ly 
values are unhkely to stem solely from redshift uncertainties. 

5. CONTRIBUTION TO THE COSMIC STAR FORMATION DENSITY 

Using the 24/im SFR estimates presented here, and the 
number density of DRGs in the MUSYC fields (known to a 
factor of two; Quadri et al., in preparation), we estimate a 
star formation rate density (SFRD) for the DRGs of '^0.02 
M0yr"'Mpc"' over the redshift range 1.5 < z < 2.5, and to a 
depth of Kab = 22.9 mag. At this depth our sample is 100% 
complete, and therefore requires no completeness correction. 
A number of uncertainties affect this estimate, including those 
described in the previous sections: errors in the photometric 
redshifts, the uncertain conversion of an observed 24/im flux 
to a SFR, and the potential contamination from AGN. Also 
important are possible differences in dust properties and the 
Initial Mass Function from low to high redshift, which are 
difficult to quantify. Nevertheless, though the uncertainties in 
the individual star formation rates are high, the average star 
formation rate of the population is known to better accuracy 
since the redshift uncertainties, at least, are largely statistical. 

The total SFRD, integrated over all galaxies, provides a 
census of the star formation history of the universe, and 
the relative contributions of different galaxy populations is 
indicative of their importance in assembling stellar mass. 
In a general sense, there is now reasonable agreement be- 
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tween IR and optical SFRD estimates at z ^ 2. iReddv et alJ 
(|2005) estimate a SFRD (from X-ray measurements) of 0.10 
M0yr"'Mpc"^ over 1.4 < z < 2.6, for optical and NIR se- 
lected galaxies, with DRGs p roducing 20% of this total. Like- 
wise, Chapman et al. (2005) estimate ^ 0.1 M0yr"'Mpc"-' 
for submm-selected systems, after extending the submm lu- 
minosity function down to ^gso^im = 1 mJy, and applying a 
correction for AGN contamination. Thus, the two techniques 
converge on a similar value, and assuming the ^gso^m ^ 1 mJy 
population encompasses the bulk of the optical and NIR se- 
lected galaxies, and vice versa, this can be taken as a total 
SFRD at z 2 (but see both papers for a discussion of uncer- 
tainties and assumptions). 

Our completeness limit of Kab = 22.9 mag is roughly equal 
to the Reddv et al. (2005) DRG depth over a similar redshift 
range, and we find the same SFRD for DRGs of ~ 0.02 
Moyr"'Mpc"''; thus, the X-ray and MIR determined SFRD 
measurements agree that the DRGs produce roughly 20% of 
the total SFRD of optical and NIR selected star forming galax- 
ies. A consistent picture is seen at longer wavelengths. If 
we use the 24/im flux density of each DRG to predict its ap- 
proximate 850/im flux density, we find that our SFRD esti- 
mate is dominated by objects with ^gso^^m > ImJy, and thus 
again DRGs produce ^ 20% of the total SFRD as measured 
at submm wavelengths. In conclusion, DRGs, to a depth of 
A'ab = 22.9 mag are important contributors to the total SFRD 
at z ~ 2, and by extension, to the assembly of massive galax- 
ies. 



6. SUMMARY 

The analysis presented here indicate that DRGs constitute 
a heterogeneous population, with the majority consisting of 
luminous and dusty starburst galaxies. We find that >65% 
of the E-HDFS sample are such system s and this is in goo d 
agreement with similar Spitzer studies ( Paoovich et al."2005|, 
and rest-frame o ptical SED modeling (Forster Schreiber et al| 
2004; Labbe et al. 2005). Using the 24-fim flux we estimate 
an average SFR for the population (to Kab = 23.2) of 130±30 
M0yr"', in line with earlier estimates from optical, submm 
and X-ray studies. Normalizing the SFR to the rest-frame 
Lv yields SFR/Ly - 19 M0yr-i(Lo x 10-'°)-', but with sub- 
stantial scatter around this value. The scatter cannot be ex- 
plained solely by redshift uncertainties, and indicates a real 
difference in the individual properties DRGs. Overall, DRGs 
are important contributors to the SFRD at z ~ 2 at a level of 
0.02 M0yr"'Mpc"^ to our 100% completeness limit of Kab = 
22.9. 
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